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Role of dense matter in collective supernova neutrino transformations 
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For neutrinos streaming from a supernova (SN) core, dense matter suppresses self-induced flavor 
transformations if the electron density Ue significantly exceeds the neutrino density n^ in the con- 
version region. If Ue is comparable to n^, one finds multi-angle decoherence, whereas the standard 
self-induced transformation behavior requires that in the transformation region n^ is safely above Tie. 
This condition need not be satisfied in the early phase after SN core bounce. Our new multi-angle 
effect is a subtle consequence of neutrinos traveling on different trajectories when streaming from a 
source that is not point-like. 

PACS numbers: 14.60.Pq, 97.60.Bw 
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I. INTRODUCTION 

The neutrinos streaming from a supernova (SN) core or 
from the aecretion torus of coalescing neutron stars are so 
dense that the neutrino-neutrino interaction causes col- 
lective flavor transformations 111 

M M 51 SI Q E m, [13, m 

2m, [22, [23, [23] • -^^ the same time, the density of ordi- 
nary matter is also very large, suppressing normal flavor 
conversions unless the neutrinos encounter an MSW res- 
onance [2a, [20, [31I, [321 • On the other hand, one of the 
many surprises of self-induced flavor transformations has 
been that dense matter barely affects them. They are 
driven by an instability in flavor space that is insensitive 
to matter because it affects all neutrino and antineutrino 
modes in the same way. Therefore, it can be transformed 
away by going to a rotating frame in flavor space [3, Q 
in a sense to be explained below. Various numerical sim- 
ulations conflrm this insight. 

We here clarify, however, that the matter density can 
not be arbitrarily large before it affects collective flavor 
conversions after all. The matter term is "achromatic" 
only if we consider the time-evolution of a homogeneous 
(but not necessarily isotropic) neutrino ensemble on a 
homogeneous and isotropic matter background. If the 
matter background is not isotropic, the current-current 
nature of the neutrino-electron interaction already im- 
plies that different neutrino modes experience a different 
matter effect. 

It is more subtle that even without a current, mat- 
ter still affects different neutrino modes differently if we 
study neutrinos streaming from a source. The relevant 



*UMR 7164 (CNRS, Univorsito Paris 7, 
Paris) 



CEA, Observatoirc do 



evolution is now the flavor variation of a stationary neu- 
trino flux as a function of distance. For a spherically sym- 
metric situation, "distance from the source" is uniquely 
given by the radial coordinate r. Neutrinos reaching a 
certain r have travelled different distances on their trajec- 
tories if they were emitted with different angles relative 
to the radial direction. Therefore, at r they have accrued 
different oscillation phases even if they have the same vac- 
uum oscillation frequency and even if they have experi- 
enced the same matter background. In other words, if we 
project the ffavor evolution of different angular modes on 
the radial direction, they have different effective vacuum 
oscillation frequencies even if they have the same energy. 
The same argument applies to matter that modifies the 
oscillation frequency in the same way along each trajec- 
tory, but therefore acts differently when expressed as an 
effective oscillation frequency along the radial direction. 



The neutrino-neutrino interaction, when it is suffi- 
ciently strong, forces different modes to reach a certain 
r with the same oscillation phase. To achieve this "self- 
maintained coherence" the neutrino-neutrino term must 
overcome the phase dispersion that would otherwise oc- 
cur. Such a dispersion is caused not only by a spectrum 
of energies, but also by a matter background. 



To clarify and quantify these ideas we revisit in Sec. |TT] 
the time evolution of a homogeneous ensemble and intro- 
duce the notion of the "rotation-averaged equations of 
motion." In Sec. IIIII we discuss the geometric modifica- 
tions when we study neutrinos streaming from a source 
and quantify the multi-angle effect caused by ordinary 
matter. We conclude in Sec. IIVI with a discussion of the 
implications of our findings for realistic situations rele- 
vant for supernova neutrinos. 



II. HOMOGENEOUS ENSEMBLE 
A. Isotropic background 

To understand the role of matter in collective neutrino 
transformations we begin with the equations of motion 
(EOMs) in their simplest form, relevant for a homoge- 
neous (but not necessarily isotropic) gas of neutrinos. We 
only consider two-flavor oscillations where the most eco- 
nomical way to write the EOMs is in terms of the usual 
flavor polarization vectors Pp for each mode p and anal- 
ogous vectors Pp for the antineutrinos, 



Pp = Hp X Pp 



The "Hamiltonian" is 



Hp = tjpB + AL 



Ai / dq (Pq - Pq) (1 - V, 



(1) 



(2) 



with dq = d'^q/(27r)'^. We have used the vacuum oscil- 
lation frequency ujp = ATn'^/2E with E = |p| for rel- 
ativistic neutrinos, B is a unit vector in the mass di- 
rection in flavor space, and L is a unit vector in the 
weak-interaction direction with B • L = cos 26* and 9 
the vacuum mixing angle. The effect of a homoge- 
neous and isotropic medium is parametrized by A = 
y2GFine- — n-e+) whereas the neutrino-neutrino term 
is given hy fj, = ^/2G-pnv^- For simplicity we here as- 
sume that initially only Vf. and D^ are present with an 
excess neutrino density of rt^^ = (1 + e) np^ [13[. The 
polarization vectors are initially normalized such that 
l/dpPpI = 1 and | JdpPpj = 1 + e. For antineutri- 
nos the Hamiltonian is the same with Up —> —cUp. 

The matter term is "achromatic" in that it affects all 
modes of neutrinos and antineutrinos in the same way. 
It was first pointed out by Duan et al. [31 that therefore 
one may study the EOMs in a coordinate system that ro- 
tates around L with frequency A so that the matter term 
disappears. In the new frame the vector B rotates fast 
around L so that its transverse component averages to 
zero. Therefore, in the new frame the rotation-averaged 
Hamiltonian is 



(H) 



ujcos{2e)L + ^iiT' -F), 



(3) 



where for the moment we consider the even simpler case 
of an isotropic and monochromatic neutrino ensemble 
where the entire system is described by one polarization 
vector P for neutrinos and one P for antineutrinos. 

A dense matter background effectively projects the 
EOMs on the weak-interaction direction. In particular, 
the relevant vacuum oscillation frequency is now uj cos 26. 
For a small mixing angle, the case usually considered in 
this context, this projection effect is not important. How- 
ever, a large mixing angle would strongly modify the pro- 
jected LU. Maximal mixing where cos 20 = would pre- 
vent any collective flavor transformations, an effect that 
is easily verified in numerical examples. 



One usually assumes that the (anti)ncutrinos arc pre- 
pared in interaction eigenstates so that initially P and 
P are oriented along L. Therefore, the rotation-averaged 
EOMs alone do not lead to an evolution. However, in the 
unstable case of the inverted mass hierarchy, an infinites- 
imal disturbance is enough to excite the transformation. 
The fast-rotating transverse B component that was left 
out from the EOM is enough to trigger the evolution, but 
otherwise plays no crucial role [7|, |9| . 

If we consider a homogeneous system where ^ is a 
slowly decreasing function of time, one can find the adi- 
abatic solution of the EOM for the simple system con- 
sisting only of P and P [ll|, [l^. In vacuum, this is a 
complicated function of e and cos 20. In dense matter, 
however, we are effectively in the limit of a vanishing 
mixing angle because the initial orientation of the polar- 
ization vectors now coincides with the direction relevant 
for the rotation-averaged evolution. The original vacuum 
mixing angle only appears in the expression for the pro- 
jected oscillation frequency uj cos 26. 

With z = Pz the adiabatic connection between P^ and 
fi is now given by the inverse function of 

w cos 26* e + 2z e + 2z + (3e + 2z)z 

= ^ ' — (4) 

M 2 2^{l + z)il + z + 2e)' 

where — 1 < z < -f 1. The "synchronization radius" rgync 
where the adiabatic curve begins its decrease is implied 
by P2 = z = 1 . One finds the familiar result 0, [ll[ ^ 



; cos 20 
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sync 



(5) 



For fi values larger than this limit, the polarization vec- 
tors are stuck to the L direction. 

Without matter one finds that (P — P) • B is con- 
served [9|. Here, the analogous conservation law applies 
to (P — P) • L. Therefore, the adiabatic solution for P^ 
is such that P^ — Pz is conserved, i.e., Pz = Pz + £■ 

In summary, the presence of dense matter simplifies 
the EOMs and in that the adiabatic solution is the one 
for a vanishing vacuum mixing angle, provided one uses 
the projected vacuum oscillation frequency. 



B. Background flux 

As a next example we still consider a homogeneous sys- 
tem, but now allow for a net flux of the background mat- 
ter, assuming axial symmetry around the direction de- 
fined by the flux. For simplicity we consider a monochro- 
matic ensemble with a single vacuum oscillation fre- 
quency LU. We characterize the angular neutrino modes 



^ In these papers fi was normalized to the density of neutrinos and 
the results were expressed in terms of a = 1/(1 -|- e). Here we 
have normalized /j, to the density of antineutrinos and use the 
picture of an excess of neutrinos, expressed by e. For the same 
physical system our fj, is the one of Ref. Q], divided by 1 -I- e. 



by their velocity component v along the matter flux di- 
rection. The EOMs are in this case 



ir y -tl^ X STy . 

The Hamiltonian for the mode v is 



(6) 



H.„ =LuB + {\- X'v)L + fi{T> - vF) , (7) 

where A' = Xve , Ve being the net electron velocity. More- 
over, 



D 
F 



dV {-Py ~Py) , 



dV V (Py — Py 



(8) 



are the net neutrino density and flux polarization vectors. 
The normalization is / du Pt, = 1 and J dvPy ^ I + e. 

Next, we transform the EOMs to a frame rotating with 
frequency A, allowing us to remove the matter term, but 
not the matter flux. 



(H„) = {u;- \'v)L + ^(D - vF) . 



(9) 



We assume 9 to be small and thus use www cos 29. We 
now have a system where the effective vacuum oscilla- 
tion frequencies for neutrinos are uniformly distributed 
between w ± A' and for antineutrinos between —to ± A'. 
Even after removing the average common precession of all 
modes, their evolution is still dominated by the matter- 
flux term if A' 3> H- In other words, collective behavior 
now requires /x > A' and not only ji^ uj. 

The simplest example is the "flavor pendulum" where 
for e = and an isotropic neutrino gas one obtains the 
well-known pendular motions of the polarization vectors. 
Matter docs not disturb this behavior, except that it 
takes logarithmically longer for the motion to start Q. 
However, a matter flux, if sufhciently strong, suppresses 
this motion and the polarization vectors remain stuck to 
the L direction for both mass hierarchies. If the neutrino 
distribution is not isotropic, the ensemble quickly deco- 
heres kinematically [i2| , an effect that is also suppressed 
by a sufficiently strong matter flux. 

We have verified these predictions in several numeri- 
cal examples, but have not explored systematically the 
transition between a "weak" and a "strong" matter flux 
because a homogeneous ensemble only serves as a concep- 
tual example where matter can have a strong influence 
on self-induced transformations. 



The neutrino density decreases with distance so that one 
can mimic this situation by a homogeneous system evolv- 
ing in time with a decreasing density, the expanding 
universe being a realistic example. However, the anal- 
ogy has important limitations because collective oscilla- 
tions introduce geometric complications into the spatial- 
variation case. 

The simplest nontrivial example is a perfectly spheri- 
cal source ( "SN core" ) that emits neutrinos and antineu- 
trinos like a blackbody surface into space. The matter 
background is also taken to be perfectly spherically sym- 
metric, but of course varies with radius. As a further 
simplification we consider monochromatic neutrinos and 
antineutrinos that are all emitted with the same energy. 
Therefore, the only variable necessary to classify the neu- 
trino modes is their angle of emission. The most useful 
variable is u = sin da where -dn is the angle of emission 
at the neutrino sphere at r = i? [ij] . At a distance r the 
radial velocity of a mode u is 



Vu,r = a/1 - u{R/rY 



(10) 



For blackbody-like emission the flux modes are uniformly 
distributed in the interval < u < 1. 

The polarization vectors are taken to represent the 
neutrino flux density because it is the flux integrated over 
a sphere of radius r, not the density, that is conserved 
as a function of r. Ignoring a possible matter flux, the 
EOMs are [11 



B P 

^r-^ u.r 

where the Hamiltonian is 
LuB + A,.L 

■n-7/,.r 



^u,r X ±u,r J 



Mr 



D, 



(11) 



(12) 



For antineutrinos we have, as always, uj -^ —lu. Since 
the polarization vectors describe the fluxes, the global 
density and flux polarization vectors are 



F, 



du ■ 



du (P„ 



P ) 

^ u,rj 1 



(13) 



using the normalization J„ du P u.r = 1 and J„ du P^.^ ~ 

1+e. The matter coefficient Ar = \/2 Gpi^e- ('')~"-e+ (^)] 
encodes the effective electron density at radius r whereas 



U' 



[i-r = [i-R- 



(14) 



III. SPHERICAL STREAM 

The most general case of neutrino flavor evolution con- 
sists of an ensemble evolving both in space and time. In 
practice, however, one usually considers quasi-stationary 
situations where one asks for the spatial flavor variation 
of a stationary neutrino flux streaming from a source. 



Here, [ir = \f2GY^ij^{K) with <i>p^(i?) the antineutrino 
flux at the neutrino sphere. Therefore, \ir always varies 
as r^^ due to the geometric flux dilution, whereas Ar is 
given by the detailed matter profile of a SN model [30, 

in, 111. 

If we ignore for the moment the neutrino-neutrino term 
in Eq. P^ . the EOM written in terms of the radial 



coordinate r simply re-paramctrizcs the distance along 
the trajectory of a neutrino. Naturally, the oscillation 
pattern of a trajectory that is tilted relative to the ra- 
dial direction produces a compressed oscillation pattern 
when projected on the radial direction. In the absence 
of neutrino-neutrino interactions, this simply means that 
along the radial direction the flavor variation of the global 
neutrino stream decoheres kincmatically, even if the neu- 
trinos are monochromatic and thus have the same vac- 
uum oscillation frequency along their trajectories. 

In the presence of neutrino-neutrino interactions, kine- 
matical decoherence among different angular modes can 
still occur and in fact is a self-accelerating process if the 
asymmetry e is too small. On the other hand, for a suf- 
ficient asymmetry, the radial variation of different an- 
gular modes is collective and they behave almost as if 
they were all emitted with the same angle relative to the 
radial direction. While this quasi single-angle behavior 
is theoretically not understood, numerically it has been 
consistently observed @, [iJ, [l3l ■ 

The variation of the polarization vectors with the com- 
mon radial coordinate r now acquires dynamical signif- 
icance in that the polarization vectors evolve differently 
than they would in the absence of neutrino-neutrino in- 
teractions. From Eq. ([72]) it is obvious that the matter 
term is no longer the same for all modes and thus can 
not be transformed away by going to a rotating frame. 
This behavior does not depend on the radial variation of 
Ar — even a homogeneous medium would show this multi- 
angle matter effect. 

For quasi single-angle behavior to occur, e must not be 
too small, a condition that is probably satisfied in a re- 
alistic SN. Therefore, the synchronization radius implied 
by Eq. ([5]) is always much larger than the neutrino-sphere 
radius i?, allowing us to expand the EOMs in powers of 
R/r < 1. Using 



where 



1'.'" ^ 2 ^2 



r2 



(15) 



we find 



H„,, = (cjB + A,L) (1,^2 



u R' 

2 72 



n2 



2r2 



where 



^r 



:,r j 



(16) 



(17) 



and Fr is the same as before. At large r a small correc- 
tion to oj is not crucial and can be ignored. The radial 
variation of A^ is slow compared to the precession, so we 
can go to a frame that rotates with a different frequency 
Ar at each radius. Finally the rotation- averaged Hamil- 
tonian is 



Mr 



M, 



R^ 

2r2 

2r2 



= y"fl 



R^ 

2r4 



(19) 



and we have assumed lj cos 29 k uj. 

The multi-angle matter effect can be neglected if in 
the collective region beyond the synchronization radius 
we have 



equivalent to 



a; < M^ 



ng+ <C np 



(20) 



(21) 



In the opposite limit we expect that the large spread of 
effective oscillation frequencies prevents collective oscilla- 
tions. In this case all polarization vectors remain pinned 
to the L direction and no flavor conversion occurs. 

For intermediate values it is not obvious what will hap- 
pen. One may expect that the multi-angle matter effect 
triggers multi-angle decoherence, destroying the quasi 
single-angle behavior. This indeed occurs for the inverted 
hierarchy whereas in the normal hierarchy we have not 
found any conditions where multi-angle decoherence was 
triggered by the multi-angle matter effect. We recall that 
for a sufficiently small e multi-angle decoherence occurs 
even in the normal hierarchy whereas no collective trans- 
formation arise for a sufficiently large e |13j . 

We illustrate these points with a numerical example 
where i? = 10 km, w = 0.3 km"\ 9 = lO'^, /x^ = 
7 X 10" km"^ and 



Ar — ^R [ — 



(22) 



(H„ 



(cj + M A*) L + M* (Qr + UF,.) 



(18) 



with n = 2. This particular value of the power-law in- 
dex leads to the same radial dependence of /x* and A*, 
see Eq. P^ . In Fig. [T] we show the radial variation /i* 
and A* for different choices of A^, between lO'^ km^^ and 
10^ km~^. Even for the smallest matter effect, the ordi- 
nary MSW resonance, defined by the condition A^ = tj, 
stays safely beyond the collective region. 

In Fig. [2] we show the corresponding variation of Pz 
for three different cases: inverted mass hierarchy and e = 
0.25 (top panel), inverted mass hierarchy and e = 0.06 
(middle panel), and normal mass hierarchy and e = 0.06 
(bottom panel). In the top panel we observe the usual 
transformation for a small matter effect, a complete sup- 
pression of transformations for a large matter effect, and 
multi-angle decoherence for intermediate cases. Repeat- 
ing the same exercise for the normal mass hierarchy and 
the same e reveals no macroscopic influence of the matter 
term. 

For a sufficiently small e one finds self-induced multi- 
angle decoherence for both hierarchies. In the middle and 
bottom panels of Fig. [2] we show how with a sufficiently 
strong matter effect the decoherence can be suppressed 
for both hierarchies. 
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FIG. 1: Radial variation of /i* and A* for our numerical exam- 
ples with the indicated value of A_r. The value of uj considered 
in this work is also shown. 



IV. DISCUSSION 

We have identified a new multi-angle effect in collec- 
tive neutrino transformations that is caused by a matter 
background. Previous numerical studies of multi-angle 
effects had used a matter profile that satisfies the con- 
dition Eq. ((2T|) in the critical region [H, [l8|. In other 
multi-angle studies matter was entirely ignored [ij] and 
otherwise, single-angle studies were performed. There- 
fore, the multi-angle matter effect discussed here had es- 
caped numerical detection. 

In many practical cases relevant for SN physics or in 
coalescing neutron stars, the density of matter is proba- 
bly small enough so that this effect can be ignored. On 
the other hand, for iron-core SNe, during the accretion 
phase the matter density can reach values and can have 
such a profile that it is important. The possibility of such 
large matter densities had led some of us to speculate 
that the second-order difference between the Vf^ and i^r 
refractive effect could sometimes play an interesting role 
for three-flavor collective transformations [2l|]. It is easy 
to show, however, that the density requirement for this 
mu-tau effect to be important implies that the multi- 
angle matter effect can not be avoided. In this sense, 
one complicated effect caused by a large matter density 
annihilates another one. 

If at early times the matter density profile is such that 
our multi-angle effect is important, this will not be the 
case at later times when the explosion has occurred and 
the matter profile contracts toward the neutron star. In 
principle, therefore, interesting time-dependent features 
in the oscillation probability can occur. 
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FIG. 2: Radial variation of P, for three different scenarios: 
inverted mass hierarchy and e = 0.25 (top panel), inverted 
mass hierarchy and e = 0.06 (middle panel), and normal mass 
hierarchy and e = 0.06 (bottom panel). In each panel differ- 
ent values of A^ have been assumed, reported in the bottom 
panel. 



A large matter effect can be "rotated away" from the 
EOMs when it is identical for all modes. Here we have 
seen that even a perfectly uniform medium provides a 
multi- angle variation of the matter effect. We note that 
the matter fiuxes would not be important, in contrast 
to our first example of a homogeneous ensemble, because 
the relevant quantity is the spread of the matter effect be- 
tween different modes. Therefore, whenever a flux term 
would be important, the matter density term already pro- 
vides a strong multi- angle effect. 



Possible polarizations of the electron-positron back- 
ground in the strong SN magnetic fields would produce 
an additional axial term in the matter potential, propor- 
tional to the scalar product of the neutrino momentum 
with the electron polarization direction [33|. However, 
we expect that also this effect should be negligible with 
respect to the matter density term. 

In addition, the medium can have density variations 
caused by convection and turbulence [3J] that is known 
to affect the MSW resonance under certain circum- 
stances [33, [3a, [33, S, |33- Density variations in the 
transverse direction to the neutrino stream lines may well 
cause important variations of the matter effect between 
different modes. It remains to be investigated in which 
way collective flavor transformations are affected. 
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